New Limits on Coupling of Fundamental Constants to Gravity 
Using ^"^Sr Optical Lattice Clocks 
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The ^Sq-'^Po clock transition frequency usr in neutral *^Sr has been measured relative to the Cs 
standard by three independent laboratories in Boulder, Paris, and Tokyo over the last three years. 
The agreement on the 1 x 10~^^ level makes usr the best agreed-upon optical atomic frequency. We 
combine periodic variations in the Sr clock frequency with Hg and H-maser data to test Local 
Position Invariance by obtaining the strongest limits to date on gravitational-coupling coefficients 
for the fine-structure constant a, electron-proton mass ratio and light quark mass. Furthermore, 
after ^^''Hg+, "^Yb+ and H, we add *^Sr as the fourth optical atomic clock species to enhance 
constraints on yearly drifts of a and /i. 
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Frequency is the physical quantity that has been mea- 
sured with the highest accuracy. While the second is 
still defined in terms of the radio-frequency hyperfine 
transition of ^'^^Cs, the higher precision and lower sys- 
tematic uncertainty achieved in recent years with op- 
tical frequency standards promises tests of fundamen- 
tal physics concepts with increased resolution. For ex- 
ample, some cosmological models imply that fundamen- 
tal constants and thus atomic frequencies had differ- 
ent values in the early universe, suggesting that they 
might still be changing. Records of atomic clock fre- 
quencies measured against the Cs standard can be ana- 
lyzed P, 0| to obtain upper limits on present-day varia- 
tions of fundamental constants such as the fine-structure 
constant a = / iAjTenhc) or the electron-proton mass 
ratio /i — rrie/rnp [1, S S B • Some unification theo- 
ries imply violation of Local Position Invariance by pre- 
dicting coupling of these constants to the ambient grav- 
itational field. Such a dependence could be tested with 
a deep-space clock mission Q, but would also be ob- 
servable in the frequency record of earth-bound clocks as 
Earth's elliptic orbit takes the clock through a varying 
solar gravitational potential (oj . Annual changes in clock 
frequencies can thus constraiii gravitational coupling of 
fundamental constants 0, [13, [13] ■ Good constraints ob- 
tained from such analyses require high confidence in the 
data and a fast sampling rate. However, a full evalua- 
tion of an atomic clock system takes several days so that 



high-accuracy frequency data is naturally sparse. 

Three laboratories have measured the doubly for- 
bidden ^^Sr ^So-"^Po intercombination line at i^sr = 
429 228 004 229 874 Hz with high accuracy over the last 
three years. These independent laboratories in Boulder 
(USA), Paris (France), and Tokyo (Japan) agree at the 
level of 1.7 Hz 12, flj, [3, El- The agreement between 
Boulder and Paris is 1 x 10~^^ [l^. ll3L[l3|. approaching 
the Cs limit, which speaks for the Sr lattice clock system 
as a candidate for future redefinition of the SI second and 
makes I'sr the best agreed-upon optical clock frequency. 
In this paper, we analyze the international Sr frequency 
record for long-term variations and combine our results 
with data from other atomic clock species to obtain the 
strongest limits to date on coupling of fundamental con- 
stants to gravity. In addition, our data contributes a 
high-accuracy measurement of an optical atomic clock 
species, which itself has low sensitivity to variation in 
fundamental constants, to the search for drifts of funda- 
mental constants, improving confidence in the null result 
at the current level of accuracy. 

In a strontium lattice clock, neutral fermionic *^Sr 
atoms are trapped at the anti-nodes of a vertical one- 
dimensional optical lattice at the Stark-cancellation 
wavelength, creating an ensemble of nearly identical 
quantum absorbers at /xK temperatures. The ^Sq- 
■^Pq clock transition [3] is interrogated with a highly 
frequency-stabilized 698 nm spectroscopy laser in the re- 
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FIG. 1: (color online), (a) Spectrum of the ^''Sr ^So-^Po clock 
transition with quality factor 2 x 10^*. (b) Measurements of 
clock transition from JILA (red circle), SYRTE (green trian- 
gle), and U. Tokyo (blue square) over the last 3 years. Fre- 
quency data is shown relative to i^o = 429 228 004 229 800 Hz. 
Weighted linear (dotted line) and sinusoidal (solid line) fits 
determine a yearly drift rate and an amplitude of annual vari- 
ation, (c) Zoom into the four most recent measurements, 
showing agreement within 1.7 Hz and giving the dominant 
contribution to both drift and annual variation. 



solved sideband limit and the Lamb-Dicke regime [IJ, 
13, 14, 15,[l3, Using individual magnetic sublevels, 
spectra with quality factors of > 2 x 10^^ have been re- 
covered [l^ as shown in Fig. [TJa). This high-resolution 
spectroscopy afforded by the optical lattice allows mea- 
surement of the clock frequency with high accuracy and 
evaluation of systematic uncertainties at one part in 10"'^^, 
limited by blackbody and residual density effects [20| . 
Spectroscopic information from the atomic sample is used 
to steer the laser to match the clock transition frequency, 
which is then measured relative to the Cs standard using 
an octave-spanning optical frequency comb 21 1 . 

In combination with data from other optical atomic 
clock species, variations in the measured Sr clock fre- 
quency can constrain variation of fundamental constants. 
It is necessary to analyze a diverse selection of atomic 
species to rule out species-dependent systematic effects 
and test the broad predictions of the underlying rela- 
tivistic theory. We will introduce the formalism required 
to constrain the coupling to gravity by first analyzing the 
global frequency record for linear drifts in a and fi. 

Figure [TJb) displays Sr clock frequency measurements 
since 2005. The frequency uncertainties are based on val- 
ues from references [H 111, Q E El, E, El, [13 . The 
date error bar indicates the time interval over which each 
measurement took place. A weighted linear fit (dotted 
line) results in a frequency drift of (— 1.0±1.8) x 10~^^/yr, 
mostly determined by the difference between the last 
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FIG. 2: (color online). The upper panel shows fractional fre- 
quency drifts for ^"Yb+ (orange), ^''Sr (red), H 1S-2S (blue) 
and ^^^Hg^ (purple) versus their sensitivity to a-variation 
relative to Cs. Sensitivity due to Cs is indicated as a dotted 
vertical line. A linear fit (solid line) determines yearly drift 
rates Sa/a and Sfi/jj,. The drift rate constraints from each 
species are shown in the lower panel as respectively colored 
bars. The fit determines a confidence ellipse (white) [3, S 0] 
with projections equal to the parameters' 1-a uncertainties. 



three high-accuracy measurements [ij, ll3|, [IJ]. This 
yearly drift can be related to a drift of fundamental con- 
stants via relativistic sensitivity constants -f^rei- Values 
for various clock transitions of interest have been calcu- 



lated in references [2J, |25| and the fractional frequency 
variation of an optical transition can be written as 
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The Cs standard operates on a hyperfine transition, 
which is also sensitive to variations in fi. For a hyperfine 
transition, the above equation is modified to 
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Here, the change in /i arises from variations in the nuclear 
magnetic moment of the Cs atom [25]. The following drift 
analysis will focus on optical clocks measured against Cs, 
since inclusion of hyperfine clock data from Rb/Cs [1,[2^ 
does not change the results significantly. 

The overall fractional frequency variation Xj of an op- 
tical clock species j compared to Cs can be related to 
variation of a and as 
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FIG. 3: (color online). Earth (blue) orbiting around Sun (red, 
mass m©) in gravitational potential U on an orbit with semi- 
major axis a, eccentricity e (exaggerated to show geometry) 
and angular velocity fi. Earth is shown at radial distance r 
from the Sun. The eccentric anomaly E is the angle between 
the major axis and the orthogonal projection of Earth's posi- 
tion onto a circle with radius a. 



For *^^Sr in particular, -cf'' = 0.06 - 0.83 - 2 = 
-2.77 [2J]. The ^'^Sr sensitivity is about 50 times lower 
than that of Cs, so that our measurements are a clean 
test of the Cs frequency variation. This allows Sr clocks 
to serve a similar role as H in removing the Cs contribu- 
tion from other optical clock experiments or to act as an 
anchor in direct optical comparisons Q . 

Other optical clock species with different sensitivity 
constants have also been analyzed for frequency drifts. 
Each species becomes susceptible to variations in both 
a and ^ by referencing to Cs. Figure [2] shows current 
optical frequency drift rates from Sr, Hg+ 7j, Yb"*" Q, 
and H Linear regression [27] limits drift rates to 



(4) 



5a/ a = (-3.3 ± 3.0) x lO^^Vyr 
5^1/ = (1.6 ± 1.7) X lO^^Vyr, 

decreasing the H-Yb+-Hg+ [1, H, errorbars (28| by 
^--^15% and confirming the null result at the current level 
of accuracy by adding high-accuracy data from a very 
insensitive species such as Sr to Fig. [2j We note that 
another limit on 5a/ a independent of other fundamental 
constants (using microwave transitions in atomic Dy ) has 
recently been reported as (— 2.7± 2.6) 10^^^/yr [29|. 

We will now generalize the formalism used for the anal- 
ysis of linear drifts to constrain coupling to the gravita- 
tional potential U and search for periodic variations in 
the global frequency record. The dominant contribution 
to changes in the ambient gravitational potential is due 
to the ellipticity of Earth's orbit around the Sun. Sup- 
pose that the variation of a fundamental constant 77 is 
related to the change in gravitational potential via a di- 



mensionless coupling constant fc^ [9]: 
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where AU{t) = U{t) — Uq is the variation in the gravita- 
tional potential versus the mean solar potential on earth 
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FIG. 4: (color online). A fit to linear constraints on gravi- 
tational coupling constants ka, and fc, from three species 
determines a plane. Its value at df^—Q, dq—Q is k^; its gradient 
along the d^ (dq) axis is k^ (kq)- The table shows sensitivity 
constants and constraints for ®^Sr, ^^^Hg"^ and the H-maser. 



Uq, and c is the speed of light. 

The variation in solar gravitational potential can then 
be estimated from Earth's equations of motion (see 
Fig. [3]). Since Earth's orbit is nearly circular, we expand 
the solar gravitational potential U{t) = —GmQ/r{t), 
with gravitational constant G, Sun mass tuq, and ra- 
dial distance Earth-Sun r(t), in the orbit's ellipticity 
e ~ 0.0167. Kepler's equation [30| relates the eccentric 
anomaly E = arccos[(l — r/a)/e\ (with semi-major axis 
a ~ 1 au) to the orbit's elapsed phase since perihelion: 



nt = E -e sin E, 



(6) 



where ~ GrriQ / ~ 2 x 10^^ s^^ is Earth's angular 
velocity from Kepler's third law. Kepler's equation has 
a solution given by a power series in the ellipticity as 
E = 57^ + 0(6), which can be used to expand 1/r and 
thus AJ7 to first order in e: 
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with a dimensionless peak-to-peak amplitude u = 
2GmQ/{ac^) ~ 3.3 x 10-i°. Thus, the ^^Sr fractional 
frequency variation due to gravitational coupling is 
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(8) 



with amplitude containing ka and fc^j as the only free pa- 
rameters. Fitting Eqn. [8] to the combined Sr frequency 
record in Fig. Wijo) gives an annual variation with am- 
plitude ysr — (—1.9 ± 3.0) X 10~^^, which constrains 
2.77ka + by division through u. 

Other atomic clock species that have been tested 
for gravitational coupling are ^^^Hg"*" Q and the H- 



maser 



ll|. H-masers are also sensitive to variations in 



the light quark mass (25|, adding a third coupling con- 



stant kq. Although the maser operates on a hyperfine 



4 



transition, the H atom is well understood, permitting 
the use of H-maser data with optical clocks to constrain 
kg . Using sensitivity coefficients from references [3, [1^ , 
each atomic clock species j contributes a constraint of 
the general form ^ 

(9) 



[ka + Cj'i^kfj + C^kq — Uj/t 



Division by gives this equation the form of a linear 
function in two variables = c^/c| and = c^/cl. 

In Fig. m each species' constraint is interpreted as a 
measurement of this linear function in the numerical co- 
efficients 



32| 



A linear fit gives: 
ka = (2.5 ±3.1) X 10"^ 
k^ = (-1.3 ± 1.7) X 10"^ 
kg ^ (-1.9 ±2.7) X 10"^ 



(10) 



Due to the orthogonal dependence on kg, the maser data 
only pivots the plane in Fig. [4] around the Hg+-Sr line, 
but its value and error bar influence neither the value nor 
the error bar of ka and fc^. The values agree well with 
zero and we conclude that there is no coupling of a, /i 
and the light quark mass to the gravitational potential at 
the current level of accuracy. We note that the coupling 
constant ka has recently been measured independe ntly 
in atomic Dy, resulting in ka = (-8.7 ± 6.6) x 10~^ 
limited by systematic effects. While optical clocks are not 
as sensitive to variations in constants as Dy, systematic 
effects have been characterized at much higher levels (2^ . 

The unprecedented level of agreement between three 
international labs on an optical clock frequency allowed 
precise analysis of the Sr clock data for long-term fre- 
quency variations. We have presented the best limits to 
date on coupling of fundamental constants to the gravita- 
tional potential. In addition, by adding a high-accuracy 
measurement of a low-sensitivity species to the analy- 
sis of drifts of fundamental constants, we have increased 
confidence in the zero drift result for the modern epoch. 
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